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ABSTRACT 

Microwave remote sensing of rough surfaces (both land and ocean), using moving 
platforms (aircraft and satellite) , as well as ground based measurements has illus- 
trated the need for a better understanding of the interaction of the radar signals 
with these surfaces. This interaction is particularly important for the ocean surface 
where the radar modulation can yield information about the long ocean wave field. 

Radar modulation measurements from fixed platforms have been made in wavetanks and 
the open oceans. The surfaces have been described in terms of two-scale models. The 
radar modulation is considered to be principally due to: (1) geometrical tilt due to 
the slope of the long ocean waves and (2) the straining of the short waves (by 
hydrodynamic interaction). For application to moving platforms, Synthetic Aperture 
Radar (SAR) and Side Looking Airborne Radar (SLAR) , this modulation needs to be 
described in terms of a general geometry for both like- and cross-polarization since 
the long ocean waves, in general, travel in arbitrary directions. In the present 
work, the finite resolution of the radar is considered for tilt modulation with 
hydrodynamic effects neglected. 

1 , INTRODUCTION 

The full wave approach is used to determine the modulation of the like- and 
cross-polarized scattering cross sections for composite models of rough surfaces 
illuminated by SAR. The full wave approach accounts for both specular point scatter- 
ing and Bragg scattering in a self-consistent manner. Thus, the total scattering 
cross section is expressed as a weighted sum of two cross sections (Bahar et al. , 
1983). The first is the scattering cross section associated with the filtered surface 
consisting of the large-scale specular components of the illuminated rough surface 
area. The second is the cross section associated with the surface consisting of the 
small-scali spectral components that ride on the filtered surface. 

Full wave solutions are derived for the scattering cross sections of a relatively 
small area or resolution cell of the rough surface that is effectively illuminated by 
SAR. The normal to an arbitrarily oriented mean plane associated with the illuminated 
cell is characterized by tilt angles ft and t in and perpendicular to a fixed reference 
plane of incidence. It is assumed that the lateral dimension of the resolution cell 
L g is much larger than both the electromagnetic wavelength and the surface height 
correlation distance for the cell. As the SAR scans different portions of the rough 
surface S, the direction of the unit vector normal to the cell F fluctuates. In this 
paper the Modulations” of scattering cross sections are determined as the tilt 
angles ft and t fluctuate . In a recent study of ’’tilt modulation” by Alpers et al. 
(1981) , first-order Bragg scatter due to capillary waves on a tilted plane is ^jnsid- 
ered. It can be shown that if the large scale spectral components of the surface 
within the -ell are ignored, the full wave solutions derived here for tilt modulation 
reduce to the results obtained by Alpers et al. 
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The expression for the physical optics (specular point) cross section for the 
large-scale surface h, is 


P 2 (rl j n 1 j n) p (n) ! 


rough surface h(x,z) and the polarization of the incident and scattered waves 
(Bahar, 1981a, b). The shadow function P 2 is the probability that a point on the 
rough surface is both illuminated and visible, given the slopes ( h^ , h ), at the 
point (Sruth, 1967; Sancer, 1969). The probability density function for the slopes 
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h x and h z is p(n). The factor x ( v ) that multiplies < 0 ^ > accounts for the degrada- 
tion of the contributions from the specular points due to the superimposed small 
scale rough surface h . 
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in which <h > is the mean square of the surface height h and v-,v~ and v_ are the 
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components of v (6) in the local coordinate system_(at each point on the large 

scale surface) associated with the unit vectors anc * Thus v can also be 

expressed as 
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The function W (v_, v__) /2^ m is the two-dimensional Fourier transform of (<h h ? >) m , 
m x 2 pQ s s 

For 6 << 1 and arbitrary p(h ,h ) the first term in (9), <0 n is also in 

X Z Si 

agreement with Valenzuela’s solutions that are ’’mostly based on physical considera- 
tions” (Valenzuela, 1968, Valenzuela et al., 1971). For small slopes n * a and 
6 << 1 the first term in (3) reduces to Brown’s solution (1978) based on a^combina- 
tion of physical optics and perturbation theory. Since it is assumed (on deriving 
(3) from the full wave solutions for the scattered fields) that the surface h* 
satisfies the radii of curvature criteria a| w§ll as the condition for deep phase 

modulation, it is necessary to choose 8 * 4k <h^> > 1 in order to assure that the 

os — 

weighted sum of cross sections (3) remains insensitive to variations in kj, the 
wavenumber where spectral splitting is assumed to occur (Bahar et al., 1983). 

In order to apply the full wave approach to SAR it is necessary to modify the 
results presented in this section (a) to account for the filtering of the very 
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For the illustrated examples presented, the scattering cross sections and their 
derivatives with respect to the tilt angles are evaluated for all angles of inci- 
dence. The modulation of the like cross sections near normal incidence is due pri- 
marily to fluctuations in specular point scattering while the modulation of the like 
cross section for near gracing angles is due primarily to fluctuations in Bragg 
scattering. Thus, for large angles of incidence the cross sections for the hori- 
zontally polarized waves are shown to be more strongly modulated than the cross 
sections for vertically polarized waves. The relative modulations of the like polar- 
ized backscatter cross sections are optimum for incident angles between 10° and 15° 
depending upon the lateral dimension of the resolution cell and the polarization. 

2. F0k: ,,, LATI0N OF THE PROBLEM 

The full wave solutions for the. normaliz ed cro?« sections per unit area are 
summarized here for composite rough surfaces. The position vector to a point on 
the rough surface is expressed as follows: 


(♦> 


r s “ r £ (x ’V z) + nh s (!) 

in which y=h^(x,z) is the filtered surface consisting of the large scale spectral 
components of the rough surface_and h g , the small scale surface height is measured 
in the direction of the normal (n) to the large scale surface y*h£. For a homogenous, 
isotropic surface height the spectral density function is the Fourier transform of 
the surface height autocorrelation function <h(x,z) ,h f (x 1 ,z 1 )> . 
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The surface h^(x,z) consists of the spectral components k - (v x + v z ) < k d and the 
remainder term h g (x,z) consists of the spectral components k > k d . The full wave 
approach accounts for both specular point scattering and Bragg scattering in a self- 
consistent manner the total scattering cross section can be expressed as a weighted 

PQ 

sum of the cross section <0 > for the filtered surface h 0 and the cross section 

PO £ * 

<o x > 3 for the surface h g that rides on the large-scale surface h (Bahar et al., 

1983 ^ <a P( ^> * < <J P( ^ > £ + < ° P< ^ > 8 • (3) 

The symbol < > denotes statistical average. The first superscript P corresponds to 
the polarization of the scattered wave while the second superscript C rresponds to 

the polarization of the incident wave. To derive (3) using the full \ approach it 

is implicitly assumed that the large scale surface meets the radii of curvature 
criteria (associated with the Klrchhoff approximations for the surface fields) as 
well as the condition for deep phase modulation. Thus the first term in (3) is 

«? PQ > jl -|x 8 (v • n g )| 2 <o PQ > (4) 

in which x 8 is the characteristic function for the small scale surface 

. X 8 (v * n ) - x 8 (v) - <exp ivh > (5) 

a “ a - -f J _ S 

v-k-k ■ k Q (n - n ) , v ■ |v| . (6) 

The unit vectors n^ and n* are in the directions of the scattered and incident wave 
normals respectively; thus for backscatter n ■ -n*. The free space radio wavenumber 
is k Q . An exp(iut) time dependence is assumed. The vector n g is the value of the unit 
vector n normal to the surface h(x,z) at the specular points. Thus 
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large scale spectral component of the rough surface by the SAR that effectively illu- 
minates a relatively small area of cell F of the rough surface S and (b) to account 
for the normal to a reference plane associated with the illuminated cell which is 
characterized by arbitrary tilt angles ft and t in and perpendicular to the reference 
plane of incidence. It is assumed e that the lateral dimension of the cell illu- 
minated by the SAR is much larger than the surface height correlation distance for 
the cell and that as the SAR scans different portions of the rough surface S the 
direction of the unit vector normal to the cell F fluctuates. Our purpose is to 
determine the Modulation” of the backscatter cross sections <c^Q> (3) as the tixt 
angles (of the normal to the cell) in and perpendicular to the reference plane of 
incidence fluctuate. 

3. SCATTERING CROSS SECTIONS FOR ARBITRARILY ORIENTED RESOLUTION CELLS OF THE 
ROUGH SURFACE 


Let x,y,z be the reference coordinate sy ^em associated with the surface of the 
cell F that is UlurMnated by the SAR such that the mean surface of the cell is the 
y=Q plane. Furthermore, let x f ,y f ,x ? be the fixed coordinate system associated with 
the large surface S such that the unit vector 3* is normal to the mean rough surface 
height hCx’jZ 1 ). The unit vector is expressed in terms of the unit vectors 

of the fixed coordinate system (x'^y^z 1 ): 

n* = -n^ = sin0 1 a 1 - cosO* a* . (13) 

ox o y 

The unit vector a normal to the reference surface associated with the cell is ex- 

y 

pressed in terms of the tilt angles ft and t in and perpendicular to the fixed 
plane of incidence, the x^y’ plane. Thus 


a = sinft cost a 1 + cosft cost a* + sint a 


(14) 


For convenience a x and a z , the unit vectors associated with the cell, can be chosen 

such that the plane of incidence in the x,y,z coordinate system is normal to the 

vector a . Thus 
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and the expression for n in the x,y,z coordinate system is 
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The angle between the plane of incidence in the fixed coordinate system (x' ,y\z') 

and the plane of incidence in the coordinate system (x,y,z) associated with the cell 
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For backscatter n*' ■ -nL Thus the angle tpp between the plane ov scatter in the fixed 
coordinate system (x',y',z') and the plane of scatter in the coordinate system asso- 
ciated with the cell is 






(18) 


The matrix that t nsforms the incident vertically and horizontally polarized waves 
in the fixed coordinate system to vertically and horizontally polarized waves in the 
cell coordinate system is therefore (Bahar, 1981a, b) 


152 


(19) 


T 


i 

F 


COSll'p 

sirup 

-sinijip 

cosip. 


Similarly, the matrix that transforms the scattered vertically and horizontally 
polarized waves in the cell coordinate system back into the vertically and hori- 
zontally polarized waves in the fixed coordinate system is 



Thus in view of (18) s T^ 
matrix D given by 
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The coefficients D 
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in (8) are elements of a 2 x 2 

( 21 ) 


in which C* n is the cosine of the angle becween the incident wave normal n^ and the 
unit vector n normal to the rough surface of the cell hp(x,z). Thus 

C in - -r^-n - cos6 in (22) 

o o 

where n* is given by (16) and n is given by (7a) with fp(x,y) * y-hy(x,y). The ele- 
ments of the scattering matrix F in (21) are functions of the unit vectors n*,n and 
n a3 well as the media of propagation above and below the rough surface S (Bahar, 
1981a). The matrix T* transforms the vertically and horizontally polarized waves in 
the cell coordinate system (3^,5 to vertically and horizontally polarized waves 

in the local coordinate system tKat conforms with the rough surface, n i» n 2» n 3 (12). 
Similarly, the matrix T^ transforms the vertically and horizontally polarized waves 
in the local coordinate system back into vertically and horizontally polarized waves 
in the cell coordinate system (Bahar, 1981a). 


To account for the arbitrary orientation of the cell, the matrix D in (21) must 
be post-multiplied by t|» and pre-multiplied by Tf . Thus the elements of the matrix 
D in (8) must be replaced by the elements of the matrix D where 


= T f D T± 


(23) 


Furthermore, in view of the effective filtering by the SAR of the very large scale 
spectral components of the rough surface f(x',z') = 0, the spectral density function 
for the rough surface fp(x,y) * 0 associated with the resolution cell F is given by 


W(v-,v_) * { 
F x z 1 


W(v x’ v z> 


k > k 
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k < k 
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where W(v_,v_) is the spectral density function for the surface S, f(x',z') * 0. 

The wavenumber k is 

k * 2 tt/L < k, (25) 

s s d 


where L s is the width of tne area of the cell illuminated by the SAR. The very large 
scale surface consisting of the spectral components 0 < k < k 8 are responsible for 
tilting the resolution cell with respect to the mean sea surface. 

Thus on replacing the spectral density function U (2a) for th** surface S by the 
spectral density function Wp for the cell F (24) and on replacing the elements 
of the matrix D by the elements D^Q of the matrix (23) the expression (3) can be 
used to determine the normalized backscatter cross section for an arbitrarily 
oriented cell F. In view of (19) and (20) the expressions for these backscatter 
cross sections are explicit functions of the tilt angles 0 and -t. For the special 
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case t * 0 (tilt is in the plane of incidence) the matrices T* and reduce to 
identity matrices and 


Thus for t - 0 
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Therefore to obtain 3<a >/3ft for Q * 0 and x ■ 0 it is sufficient to evaluate <o x > 

as a function of 0* with both 0 and x set equal to zero. The value for 3<o^>/3x can 

either be evaluated analytically since d|Q (23) is an analytic function of x, or the 

derivative could be evaluated numerically. 

4. ILLUSTRATIVE EXAMPLES 


For the illustrative examples presented in this section, the following specific 
form of the surface height spectral density function is selected (Brown, 1978) 
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x* z ff x z 
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where W is the notation used by Rice (1951) and S is the notp.cion used by Brown (1978), 
For the assumed Isotropic model of the sea surface 
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i. which k c is the spectral cutoff wavenumbc^ (Brown 1978) and V is the surface wind 
speed. The wavelength for the elccc* cmag-T."' -v.ve is 

A » 3.0 cm (k ** 2r.,\ (cm) 1 ) . (3] 

0 0 

The relative complex dielectric coefficient for the sea is 

e r - 48 - i35 (33 

and the permeability for the sea is the same as for free space (w^ « 1) . 

The mean square height for the small scale surface h is given by 
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The mean square slope for the large scale surface h^ within the resolution cell, is 
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in which k s is given by (25) . The mean square height for the large scale surface 
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For B = 4k 2 <h 2 > = 1.0, k = 0.201. For L =300, 1000 and 2500 cm (25) 
os d s 


and 
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k Q <h^> « 21.9, 173 and 357 respectively. The slope probability 


density function within a resolution cell is assumed to be Gaussian; thus 
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and the physical optics (specular point) backscatter cross section iu (8) (Bahar, 
1981a) 
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in which 6pq is the Knnecker delta and Rp(P=V,H) is the Fresnel reflection coefficient 
for the vertically or horizontally polarized waves (Bahar, 1981a, b). 

W w W 

In Fig. la, and b <a >, and ~(d<o >/d Ji)/<o > are plotted for 0 = 0 and t = 0 as 

functions of 0^ the angle of incidence with respect to the fixed reference system 
(x'jy'jZ*). In these figures Lg ■ 300, 1000 and 2500 cm. 



VV 

Figure la. <0 >, for ft - 0 and 

T ■ 0 as a function of 6'. 

o 

(A) L fl * 300 cm, (Q) L g = 1000 cm, 

(0 ) L » 2500 cm. 

s 
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Figure lb. -(d <o >/dti)/«j > for 

n ■ 0 and t * 0 as a function of 

0 1 • ( A) L * 300 cm, 

o s 

(O) L s * 1000 cm, 

(□ ) L * 2500 cm. 
s 


HH 

In Fig. 2a, and b these results are repeated for <o >. It is interesting to note 
that the effective filtering of the very large scale spectral components of the rough 
surface (0 < k < k g ) by the SAR does not significantly change the value of unless 
L g < 300 cm. As one may expect, the modulation of the scattering cross sections in 
the plane of incidence |d<o'^>/dft| is strongest for the SAR corresponding to the 
narrowest effective beam width L g * 300 cm. Except for near-normal incidence the 
relative modulation |d<o^' > /dfl|/<a^Q> is larger for the horizontally polarized waves 
than for the vertically poxarized waves. The largest relative modulation of the like 
polarized cross sections occurs in the transition region where the contribution to 
the cross section due to Bragg scatter becomes larger than the contribution due to 
specular point scatter namely at about 10°-15° (see Figs, lb and 2b). 



HH 

Figure 2a. <o >, for ft - 

t ■ 0 as a function of O'. 

o 

(A) L - 300 cm, (O) L - 

8 S 

(□) L - 2500 cm. 
s 


0 and 
1000 cm, 
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Figure 2 b. -(<a I1H >/dQ)/<a HH > for 


^ ■ 0 and t * 0 as a function of 


e'. (A) - 300 cm, 


(O) L - 1000 cm. 


(□> L c 


5. CONCLUDING REMARKS 


, The wave a PP r °ach is used to determine the scattering cross sections for 

aroitrarily oriented resolution cells on random rough surfaces illuminated by 
synthetic aperture radars. The purpose of this analysis is to determine the modu- 
ation of the like polarized scattering cross sections as the normal to the cells 

fr 11 iY n ? P er P« n J lc i| lar t0 the plane of incidence. The full wave approach accounts 
for shadowing and both specular point scattering as well as Bragg scattering in a 
8 * 4 f “consistent manner. Thus, the scattering cross sections are expressed as 
weighted sums of two cross sections. The first cross section is associated with the 
iltered surface consisting of the large-scale spectral components of the rough 
surface. The second cross section is associated with the surface consisting of the 
small-scale spectral components. It can be shown that if the large-scale spectral 
components of the surface of the cell are neglected, the second cross section 
accounts for first order Bragg scattering and the results are in agreement with 
earlier published results (Alpers et al., 1981). However, for typical terrain or 
sea surfaces, the large-scale spectral components are not negligible. 


By using the full wave analysis, the modulation of the like and cross polarized 
cross sections can be determined for all angles of incidence and tilt angles. On 
he other hand, first order Bragg scatter theory does not account for backscattering 
near normal to the surface of the cell (Alpers et al., 1981). The results based on 
the two-scale model indicate that the relative modulation of the like polarized 
backscatter cross section is maximum for angles of incidence between 10 ° and 15° 
(depending on polarization and effective width of the resolution cell, L ). The 
analyses based on first order Bragg scatter do not provide these results! It is 
also shown that as the angle of incidence approaches zero, the modulation of the 

scattering cross sections in and perpendicular to the plane of incidence becomes 
comparable. 


When the normal to the cell is tilted in the direction normal to the plane of 
incidence (t i 0), the full wave analysis not only accounts for the change in the 
local an le of incidence 0 Q but also takes into account the fact that the local 


_ e _ _ £ . . . . U k witc * V. (.U«l It 1C iO 

planes of incidence (or scatter) are not parallel to the reference planes of 
incidence for scatter), namely ^ + o. 


Since Alpers et al. (198i) do not account for the effects of the large scale 
spectral components of the surface within the resolution cell the results presented 
here for the modulation of the like polarized scattering cross sections near 
normal incidence are significantly different from those given by Alpers et al. 
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